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The unusual GRB 110709B triggered Swift/BAT twice, with a time difference of ~ 11 minutes. Its 
light curve presented three noticeable peaks but only two were originally identihed. In this work, 
we describe each peak as due to a different central-engine phase: the hrst one is the millisecond- 
protomagnetar stage, the second one is the BH-formation collapse phase and the last one is the 
Collapsar scenario with a Blandford-Znajek engine. Additionally, we analyze and explain the 
afterglow phase evoking the standard hreball model. Our model can successfully describe the 
timescales, fluxes and spectral indices observed for GRB 110709B. 
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1. Introduction 

The Collapsar model [1,2] for long gamma-ray bursts (GRBs), requires the collapse of a 
stellar core with a specific angular momentum of ~ 10*^ ^ cm^ s^^ [3]. Under this condition the 
stellar core collapses to a black hole (BH) surrounded by an accretion disk. This implies pre¬ 
collapse spin periods of Pspin < 0.5 day (and thus, orbital periods in tidally-synchronized binaries) 
which restricts the stellar evolution prior to collapse [4, 5, 6, 7, 8, 9]. On the other hand, Duncan & 
Thompson [10] and later works [1 1] have proposed utilizing the spin down of ms protomagnetars as 
the central engines for long GRBs. These engines have much less available rotational and binding 
energy when compared to collapsars (10^^ erg vs lO^'^ erg, respectively) as the mass and spin of the 
compact object are both smaller; thus, they must be much more efficient in their energy conversion. 
Now, it is likely that magnetars may be produced by less massive stars than BHs, then they may 
be much more common. Nonetheless, it is still necessary for the progenitor of ms magnetars to 
rotate extremely rapidly prior to the collapse, both, to explain its ms rotation and, perhaps even, to 
amplify the magnetic field (e.g., [9]). In fhis work we sfudy fhe possibilify, previously suggesfed in 
Zhang el al. [12], lhal GRB 110709B was fhe resulf of a combinalion of factors lhal allowed bolh 
mechanisms, ms magnefar and collapsar, lo play a pari in fhis Iransienl evenl. 


2. Model 


To explain fhe Ihr'ee peaks (or episodes) of GRB 110709B as fhe collapse of a single slar we 
show fhe model as a funclion of lime scales, magnelic fields, masses, efc. 


First Episode; We consider a collapsing 3.5M0 core {Rpe ~ 10^ cm). The free-fall fimescale of 
such an objecl is of fhe order of 


_K I Rl^ ^ 371 / R y/V M \ ^ 

~ 2^ GMFe~ 4 VlO^cmy y.SMeJ ^ 


( 2 . 1 ) 


fhe viscous fimescale. Tv, is 

/ y \2 4Tf f 

Tv ~ ^ ~ 900s, (2.2) 

\hJ a 

and fhe rofalional kinelic energy of a profoneulron sfar (PNS; or a neufron sfar, NS) is Ihen of fhe 
order of 


Ek = 







10^cm\ g 


(2.3) 


where a Belhe, 1 B = 10^^ erg. This energy be lapped Ihr'ough a lorque exerted by fhe magnelic 
dipole. The power can be eslimaled from: 
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~ 2.2 
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3x io'^g; 
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Bs 


(2.4) 


where R is fhe radius of fhe magnelic dipole and Q. is fhe angular velocily (see Usov [13] for more 
defails). Afler a few lens of seconds, fhe ms-magnefar engine slows down and fhe jefs are shuf 
down. Mailer held oul by propeller effecl in an accretion disk will slarl slreaming down, burying 
fhe magnetic field and increasing fhe mass of fhe NS. 
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Second Episode: By assuming that around {k ~) 10% of the rest mass of the NS is released as 
binding energy during the conversion into a BH, the total energy available for the 2nd episode can 
be estimated from 

= (2.5) 

Under normal SN conditions, over 99% of the energy released leaves the star as neutrinos, without 
further interaction. For SN 1987A [21] the kinetic energy was around Ekm ~ 1 B and the total 
energy released (in neutrinos) was Ej ~ 300 B. Assuming somewhat larger efficiency (larger den¬ 
sity and temperature, thus, larger neutrino cross section), and considering we have more mass in 
the compact object we estimate the kinetic energy to be a few times larger than in SN 1987A, i.e., 
Ekin ~ 5 B, which coincides with the observations for GRB 110709B. Given the available energy 
and the fact that the material that forms the collapsing magnetar is denser than that of the accretion 
disk, a quasi-thermalized, sharp, hard-X-ray signal would be expected. 


Third Episode: After the BH forms, it accretes material and angular momentum. As the rest of 
the stellar core collapses, differential rotation and convection regenerate a strong magnetic field. 
The Kerr BH interacts with the accretion disk via the Blandford-Znajek mechanism (BZ; [14]). 
The energy can be estimated from 


Ebz 


where ai, = Jc/(GM^) = a/M, 


1,800 en/(a*) 



/(«*) = 1 



( 2 . 6 ) 


(2.7) 


and ea = 0.5, the maximum efficiency for energy extraction. The power of the BZ central engine 
is given by 

(2.8) 

(see [15, 16] for more detail). 


3. Afterglow (AG) 

In the framework of the external shocks of the fireball model the deceleration radius and 
timescale are 

^ (3.1) 

\471 mpj 

and 

' (l+z)r-^/3TJ-‘/3£l/3^ (32) 

\527lmp J 

respectively, where E is the energy, F is the bulk Lorentz factor, rj is the interstellar medium (ISM) 
density and mp is the proton mass. The AG emission is more likely to be of synchrotron origin, 
thus, the spectrum has a break in the fast cooling regime resulting in a spectral index of 1.5. 
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Figure 1: Left pane/.BAT count rates (upper panel) and photon index evolution (lower panel) during 
episodes 2 (in magenta) and 3 (in black) of GRB 110709B. The spectral model is a simple power law. 
The time is relative to the trigger time of episode 1, of 331940612 s (in MET seconds). Right panel: Fit of 
the Swift/BAT spectrum with a PL model between 835 s and 865 s after the first trigger time. The lower 
panel shows the residuals. 

4. Observations and Data Analysis 

GRB 110709B triggered the Burst Alert Telescope (BAT, [17]) onboard the Swift satellite at 
21:32:39 UT on July 9, 2011 [18]. This episode extended up to 55 s after the trigger [12]. Hereafter 
all events will be measured with reference to the first trigger time. Surprisingly, there was a second 
BAT trigger at 21:43:25 UT on July 9, 2011 after ~ 11 minutes. The emission period extended up 
to 865 s [18]; this period shows a first bump (episode 2), beginning at ~ 550 s and lasting about 60 
s (ending at 610 s), followed by a second multi-peaked bump (episode 3, from 610 s to 750 s; see 
Fig. 1) of longer duration. 

Zhang et al. [12] performed a time-resolved spectral analysis of this second emission (episodes 
2 and 3), dividing its time period in time slices of 50 s or more; they showed that the spectra 
can be fitted with cutoff power laws and that there was a strong hard-to-soft spectral evolution. 
However, their choice of the time intervals is not suitable to look at the evolution of the spectral 
parameters during episode 2, as the corresponding time slice covers its whole duration. Therefore, 
we performed a more detailed analysis of the spectra of episodes 2 and 3, by considering sub¬ 
intervals of 10 s each. We processed the S'wf/i/BAT data using the standard FTOOLS package 
(Heasoft, version 6.15). We analyzed the spectra using two different spectral models: BB-i-PL and 
PL. We found that all the spectra are well modeled with PL, while the BB-i-PL can be discarded, 
as it is not well constrained. The lower left panel in Fig. 1 shows the time evolution of the photon 
index of the PL model: it can be seen that a discontinuity in the hard-to-soft evolution comes out 
at the beginning of episode 3. This could suggest a different emission mechanism for episodes 2 


4 






















































GRB110709B 


Nissim Fraija 


and 3. We obtain that the BAT band (15-150) keV fluences of the first (from —28 to 55 s), second 
(550 to 610 s) and third (610 to 750 s) episodes are 9.54+Q}g x 10^^ erg cm^^, 2 . 31 +qq 5 x 10^^ 
erg and 8.811qj2 x 10^^ erg cm^^, respectively. 

We also analyze the spectrum between 835 s and 865 s corresponding to the last peak observed 
in the BAT light curve [12], tens of seconds after the end of episode 3. We find that the BAT 
spectrum (see right panel in Fig. 1) is well modeled with a PL with spectral photon index = 
1.59('(Qj4. By considering the value of spectral photon index, the deceleration radius (eq. 3.1 ) and 
time scale (eq. 3.2), we interpret that this peak is related to the AG of episode 3. 

5. Discussion and Conclusions 

Considering the values of the first episode: Energy E = 20 B, redshift z = 0.75 (~3.1 Gpc; see 
[20]), ISM density rj = 3 cm^^ and Eorentz factor of E = 160 and E = 50 we obtain deceleration 
radii of ~ 3.5 X 10^^ cm and R^ ~ 7.5 X 10^^ cm which correspond to time scales of G 40 s 
and G ~ 870 s, respectively. 

Taking into consideration the values for episode 3, energy E = 18.5 B, a mean ISM density T]= 
10^^ cm^^ and a bulk Eorentz factor of r= 225, we obtain a value of deceleration time G ~ 225s 
which is in agreement with the observations. 

The first jet propagates out of the star and into the ISM clearing a path, thus lowering rj from 
3 cm^^ to 10^^ cm^^ . The first AG took place at 36 s (ai = 1.55 ± 0.05 ; [12]), and the AG 
for 3rd episode occurred at 850 s (a^ = 1.59 ± 0.15) [19]. Zhang et al. (2012) [12] reported 
a component modeled with a power law of spectral index ai = 1.55 ±0.05 at the very end of the 
prompt emission, between ~ 36 and ~ 45 s which strongly suggests an early AG probably powered 
by a magnetized outflow [22, 23, 24]. 
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